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ABSTRACT 

Minkowski Functionals (MFs) are topological statistics that have become one of many standard 
tools used for investigating the statistical properties of cosmological random fields. They have 
found regular use in studies of departures from Gaussianity in a number of important cosmological 
scenarios. Important examples include the Cosmic Microwave Background (CMB), weak lensing 
studies, 21cm surveys and large scale structure (LSS). To lowest order the MFs depend on three 
generalised skewness parameters that can be shown to probe the bispectrum with differing weights. 
Recent studies have advocated the use of a power spectrum associated with the bispectrum, called 
the skew-spectrum, that has more power to distinguish between various contributions to the bis- 
pectrum than the conventional formalism adopted when using the Minkowski Functionals. In this 
article we review the motivations for studying non-Gaussianity and emphasize the importance of 
the momentum dependence of higher order correlators in investigating both inflationary and early 
Universe models as well as analytical models for gravitational instability. We then introduce the 
skew-spectra, applied to galaxy surveys, as a tool for investigating various models for primordial 
and gravitationally induced non-Gaussianities. We present analytical expressions for the skew- 
spectra for the density field and divergence of the velocity field in 3D and for projected surveys 
as a function of redshift and a smoothing angular scale. A Gaussian window function is assumed 
throughout this paper. Analytical results are derived for the case of gravitationally induced non- 
Gaussianity. These results can be generalised to incorporate redshift space effects. This will be 
useful in probing primordial and gravitationally induced non-Gaussianity from ongoing and future 
galaxy surveys. 

Key words: : Cosmology- Large Scale Structure - Minkowski Functionals - Methods: analytical, 
statistical, numerical 



1 INTRODUCTION 

The canonical model for inflation is a single free scalar field obeying the Hamilton- Jacobi slow-roll constraints in a semiclassical background 
governed by Einstein gravity with initial conditions set by assuming that the vacuum state asymptotically approaches the Bunch-Davies vacuum. 
The appeal of invoking such an inflationary model is that it provides a rather natural mechanism for solving the causality and flatness problems 
whilst allowing for the generation of primord ial density perturbations through quantum mechanical fluctuations (Mukhan ov and Chibisov|[l98ll : 
IStarobinksvll 19821 ; lHa wking 1982; Guth 1982). The scalar field model draws motivation from concepts in particle physics and in many ways rep- 
resents a vanilla model for inflation in the sense that we have made rather minimal assumptions about the presence of physics beyond the standard 
model in the early Universe. It should be emphasized, however, that whilst the scalar field models draw on ideas from particle physics there i s, as of 
yet, no known candidate particle that could source the hypothetical quantum field driving inflation (see, for example. lLvth & Riottdfl999l ), lLind"e1 
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d2005b and Mazumda r & Rochet] J20 1 lb for reviews of particle physics models for inflation). The energy scales typically associated with infla- 
tionary cosmology are currently far beyond that accessible to terrestrial experiments and we must necessarily extrapolate current theories to these 
energy scales prompting quest ions over the validity of the effective field theories in describing inflationary models and possible quantum gravity 
(i.e. U V) inspired corrections dEasfher et alj|200ll ; iMartin & BrandenbergeifeoOll ; lDanielssonll2002l ; iKachru et al]|2003l ; iBaumann & McAllistej 
l2009h . 

The CMB probes the early Universe in a linear regime acting, currently, as one of the most direct probes of the dynamics of the early Uni- 
verse. One of the largest problems with current and future CMB surveys (e.g. polarization) is the relatively small signal-to-noise ratios. Planck 
is expected to exhaust the information content of the temperature anisotropies but will not provide strong statistical constraints on the polar- 
ization with a number of dedicated polarization missions in the development stage (e.g. CMBPol, COrE). Given that we only have access to 
one realization of the Univers e, it is important to u s e as many cosmological data s ets as possible to analyse the st atistical properties of primor- 
dial density perturbations (e.g. iMunshi et alj J201 lh ; [Mead. Lewis and Kind d20lTh ; iNicholson a nd ContaldJ d2009h ). A number of recent papers 
have suggested that the statist i cs of d ark halos may pro v ide as stringent const r aints on primordial non-Gaussianity as CMB observations (e.g. 
ICarbone, Verde and Matarresd d2008h ; iDalal et alj d2008h ; iMatarrese and Verde] (2008)). When observing large scale structure we encouter more 
severe non-linearity and we find that there are a number mechanisms for generating significant non-Gaussianity. The first and most prominent is 
the non-linear evolution of structure through gravitational instability. A second source of non-linearity is the bias relation between the galaxy and 
matter distributions that arises when considering large scale structure surveys. The final, and in our case most interesting, source of non-Gaussianity 
is primordial in origin. In this paper we consider the role of large scale structure (LSS) surveys in providing a probe of early Universe cosmology 
and further develop the skew-spectra as a tool to distinguish between various contributions to the observed bispectrum of LSS. 

The MFs have been extensively developed as a statistical tool in a cosmological setting for both 2-dimensional (projected) and 3- 
dimensional (redshift) surveys. The MFs have analytically kn own results for a Gaussian random field making them suitable for studies of 



non-G a ussianity. Examples of such studies include CMB data (Natoli et al. (2010); Hikaee et al. (2008); Novikov, Schmalzing and Muk hanovl 



d2000h ; ISchmalzing & GorskH dl998h). weak tensing dMatsubara and Jainl d200l[); ISato et alj d200lh: iTaruva et al.l d2002h; IMunshi et alj d201 lh) 
large s c ale structure dGott et all dl986t);|Colesl dl988h;lGott et al.l dl989h;lMelo^dl990h ; lMoore et alJdl992h;lGott et alj dl992h;ICanavezes et al 



1998); Schmalzing & Diaferio (2000); KerscheretaU d200ll); IPark et alj d2005h ; iHikage et all d2008h; IHikage. Komatsu & Mastubar; 



2006) ; iHikage et alj d2002h). 21cm dGleser et alj d2006h ). Sunyaev-Zel'dovich (SZ) maps ( dMunshi etall 1201 lh ) and N-body simulations 
Schmalzing & Diaferid (2000); lKerscher et alj d200lh ). Note that this is an incomplete list of references and we have selected a sample of rep- 
resentative papers from the literature. The MFs are spatially defined topological statistics and, by definition, contain statistical information of all 
orders. This makes them complementary to the polyspectra methods that are defined in Fourier space. It is also possible that the two approaches 
will be sensitive to different aspects of non-Gaussianit y and systematic effects although in th e weakly non-Gaussian limit it has been shown that 
the MFs reduce to a weighted probe of the bispectrum dHikage. Komatsu & Mast ubara (200(J)). 

The next decade should see the next generation of large scale str ucture surveys beg inning to p roduce cosmological data sets in un precedented 
detail with the most notab le of these surveys including: WigglezQ ( iBlake et afcoilh . Euclid □ dLaureiis et al]|2009l) . WFIRST0 dGreen et all 
l201lh . BossQ, BigBossQ dSchlegel et alll2010h and Subaru PFS0. These future surveys can be combined with the latest CMB measurements to 
effectively trace the growth of structure formation from the surface of last scattering (z ~ 110 0) through to the local Universe ( z ~ O(l))- The 
availab ility of such data sets promises to provide tight constraints on cos mological parameters |Wang, Sper gel and Str auss ( 1999); Tegmar k et alj 
d2004ah ), angular diameter distance derived from BAO measurements dPercival et alj d2007h ). Hubble expansion rate, growth rate of structure 
forma tion dPercival and White] |2009)), Gaussianity of initial conditions and the nature of late time cosmic acceleration dDark Energy Task Force! 
(2006)). With more comprehensive data sets, a deeper understanding of systematics (e.g. Hikage et al. ( 201 1 )) and improved understanding of both 
primordial and late time contributions to non-linearity of the growth of structure (e.g. lBernardeaul d2002 ); Chen ( 2010)) it becoming possible to do 
comprehensive studies utilising higher order statistics. In particular we are slowly reaching a point where it will be possible to break degeneracies 
between various contributions to observed non-linearity (primordial, gravitational instability, bias, systematics, etc) and constrain theories for the 
early and late Universe. 



1 http://wigglez.swin.edu.au/site/ 

2 http://euclid.gsfc.nasa.gov/ 

3 http://wfirst.gsfc.nasa.gov/ 

4 http://cosmology.lbl.gov/BOSS/ 

5 http://bigboss.lbl.gov/ 

6 http://sumire.ipmu.jp/en/ 
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This paper is organised as follows. In Sj2]we briefly review the concept of Minkowski Functionals in 3D and 2D as a morphological descriptor. 
In $3]we introduce the skew-spectra. The next section [j4]is devoted to gravity induced non-Gaussianity as well as primordial non-Gaussianity. The 
section ^introduced the skew-spectra for a 3D field in general and galaxy distribution in particular. Next, in section $6]we develop the theory for 
skew-spectra for the divergence of 3D velocity field. The projected skew-spectra are considered in |7] Finally Ej9] is devoted to discussion of our 
results and conclusion. 

The particular cosmology that we will adopt for numerical study is specified by the following parameter values (to be introduced later): 
Qa = 0.741, h = 0.72, Q b = 0.044, ft C DM = 0.215, Q M = fi 6 + «cdm, n s = 0.964, w = -l,w a = 0, <r 8 = 0.803, tt v = 0. 



2 MINKOWSKI FUNCTIONALS: INTEGRAL GEOMETRIC TOOLS AS A PROBE OF NON-GAUSSIANITY 

In this section we provide a brief introduction to integral geometry and define the Minkowski Functionals for a d-dimensional space. Integral 
geometry provides a natural framework for the set of morphological descriptors of a given random field. These descriptors are intrinsically defined 
in the spatial domain and naturally take into account all n-point correlators up to arbitrary order. This provides a nice complementarity between 
the Fourier based methods and the spatially defined topological statistics. Hadwiger's characterisation theorem shows that a linear combin ation of 
a set of d+1 functionals will provide a complete morphological description of our random field in a d-dimensional space dHadwigenl 1 9591) . These 
functionals are more commonly referred to as the Minkowski functionals (MFs) and are defined, adopting an integral geometric approach, over a 
convex ring (i.e. any finite union of compact, convex sets fCd in a d-dimensional Euclidean space E d ). The importance of Minkowski Functionals 
is that they will be the unique morphological descriptors that obey motion-invariance, convex-continuity and additivity. These properties are 
desir able in a morphologica l descr iptor due to the interpretation of these properties when working on, for example, pixelised data sets (see 
e.g. Iwinitzki and Kosowksvl dl998l) ). The Minkowski functionals are commonly calculated using volume weighted curvature integrals and the 
appea l of this formalism is that there exist exact analytical results for an underlying Gaussian random field ( Adler 19811: lTomitalll986l : lGott et al.l 
1990). More recently the analytic al results have been ext e nded to weakly non-Gaussian random fields throug h a perturbative approach based 
on an Edgeworth expansion (e.g. iMatsubarj dl994l Il995h : iMatsubara & Yokoham j dl996l) : iMatsubarj {2002); Hikage, Koma tsu & Mastuban] 
(2006)) as a function of a set of skewness parameters. These results allow us to use the MFs as a test for non-Gaussianity in the perturbatively 
weak regime. This is convenient as observations from the CMB and LSS constrain non-Gaussianity to be small and, importantly, the levels of 
non-Gaussianity generated by models for inflatio n, reheating, preheating, secondary ef fects and gravitational instability are typically well within 
the domain of validity for the perturbative results (Hikage, Komatsu & Mastubara 2006). 



The MFs can be defined for both 2D (projected) and 3D (redshift) data sets. The MFs have been applied to CMB data sets dColesI 


1988; 


Park & Gott 


1988; 


Collev, Gott and Park 1996; Novikov, Schmalzing and Mukhanov 2000; Komatsu et al. 2003; Eriksen et al. 2004; Hikage et al. 


20081: Natoli 


etalJ 


20ld). large scale structure dGott et al. 1986L 19891: Melottl 199(J: Coles & Plionisl 199 ll: Moore et alj 1992; Gott et alj 1992; 


Meckeetal. 1994; Rhoads, Gott & Postman 1994; Schmalzing & Buchert 


1997; Canavezes et al. 


199 


$ Sab 


ni, Sathyaprakash and Shandarin 
ae. Komatsu & Mastubara 2006; 
aatial domain, provide a probe of 


1998; Schmalzing et al.ll 19991; ISchmalzing & Diaferio 2000; Park et al J 120051; Iffikage et alj|2002t 120031; iHika 
Hikage et al. 2008) and weak lensing surveys (Munshi et al. 201 1). The MFs, being intrinsically defined in the s 



all orders of correlation functions in contrast to the more convention al polyspectra or Fourier-space meth ods (e.g. power spectrum, bispectrum 
and trispectrum). In the weakly-non-Gaussian limit it was shown ( Hikage, Komatsu & Mastubara 2006) that the MFs reduce to a weighted 
probe of the bispectrum given in terms of a set of skewness parameters. This makes the MFs complementary as it offers an alternative probe 
of th e data in the presence of contaminants such as survey masks, inhomogeneous noise, foregrounds, etc ( Hikage, Komatsu & Mastubara 
2006). Whilst the CMB offers, perhaps, the cleanest probe of non-Gaussianity large scale structure surveys offer a complementary data set that 
can be used to investigate the nature of primordial non-Gaussianity and the nature of gravitational instability associated with the growth of structure. 

In this paper we will consider the Minkowski Functionals for 3D random fields as well as 2D projected surveys in the context of large scale 
structure. A 3-dimensional cosmological density field can be treated as a smoothed scalar field ty(x) with zero mean ('fy(x)) — 0. The excursion 
set, Q, can be calculated over a given threshold, v, by taking all the points of the random field above the threshold: Q = {x|*I'(x) > vo}, where 
a = (x]/ 2 ) 1 / 2 is the standard deviation of the field. For a three-dimensional space the Minkowski Functionals can be defined using curvature 
K2) weighted integrals: 

Voiy) = ^ [ d 3 xQ(vcr Vi(u) = / ds; V 2 (u) = f ds [«i(x) + « 2 (x)] ; Va(y) = / ds«i(x)« 2 (x). (1) 

V J v SV J dQ 6nV Jq Q AtvV Jq Q 
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The formula for MFs in both a Gaussian dTomita|["l986l) and a weakly non-Gaussian field (Mats ubarair2002h are analytically known allowing us to 
compare observational and numerical results to expected values. The generalisation to a 2D field is trivial and follows the same approach but of 
reduced dimensionality. 

Following the notation of Hikage, Ko matsu & Mastubaral d2006h . the formula for a MF in a weakly non-Gaussian field is given as a function 
of the threshold, v, with an amplitude set by Ak'. Vk{y) — AkVk{v) The amplitude of the Minkowski Functionals are weighted by the spatial 
dimension, the variance of the random field, oq, and the variance of the derivative field, a x : A k = (l/(27r) (fc+1)/2 )(cj 2 M>-fc^fc) ((Ti/y/2er ) ' 
with cjfc = 7r fc//2 /r(fc/2 + 1) being the volume of a k-dimensional unit ball. The normalised Minkowski Functionals, v k (v), can be perturbatively 
expanded with respect to the variance, <r, resulting in a Gaussian term and the corresponding non-Gaussian contributions consisting of higher order 
corrections weighted by the variance: Vk [u) = v k G ^ (y) + (v) . The Gaussian term is given by v k G ^ = e~" ^ 2 Hk-i (u) . Working to order 

<7o the non-Gaussian corrections are given by: 



V (NG) = e -^/2 



\s (0) H k+2 {v) + ~S w H k (v) + k{k 6 ^ ] S (2) H k . 2 {u) 



a + O{ol) \ , (2) 



where H k (v) are the k-th Hermite polynomials and 8^ are a set of skewness parameters. The skewness parameters are given by dMatsubari] 
l2002h : 

o-q 0o°i °~i J 2tt 2 

The Fourier coefficients are defined in terms of the density contrast of galaxies <5(x, z): <5(x, z) = l/(27r) 3 J d 3 k5(k)e lk x . Next we show how 
to generalize the skewness parameters to the corresponding skew-spectra. The skew-spectra carry more information about the relevant bispectra 
from which they are constructed. 



3 GENERALISED SKEW-SPECTRA 

Previous work has resulted in the development of optimal 3-point estimators for simple models of the bispectrum, most notably for the phenomeno- 
logical local configuration parameterised by /^f;. Most of these optimal approaches compresses the data into a single estimate based on /nl and 
consequentially we lose sensitivity in distinguishing between various contributions to the observed non-Gaussianity. 

The skew-spectra are a proposed method utilising cubic statistics constructed from the cross-correlation of two differing fields. These fields 
are constructed as a product of the maps and their derivatives. The three different skewne ss parameters that we generate will reduce to a weighted 
probe of the bispectrum. Following the method developed in iMunshi & Heavens! g2009f) we can define a power spectrum associated with each 
of these skewness parameters and hence we can associate a power spectrum to each of the MFs. The power spectrum associated with a given 
MF will, by construction, have th e same correspondence with the various skew-spectra S, as the MFs have with the one-point cumulants 
jMunshi. Smidt and Cooray||20ld) . 

The advantage to the associated power spectra is that we retain more of an ability to distinguish between various models for non-Gaussianity 
rather than collapsing all information into a single estimator (e.g. estimators for /nl). This is due to the generic momenta dependence of higher 
order correlators such as the bisp ectrum. A previous paper jMunshi. Smidt and Cooravll201(j|) derived analytic results for topological statistics 
based on the use of skew-spectra JMunshi & Heavens 2009D allowing us to relate the analytic skewness parameters to the topological properties 
of large scale structure. Each of the generalised skewness parameters or generalised cumulant correlators can be constructed from triplets of field 
variables that relate to either the original density contrast or to variables constructed from derivatives of the fields: 

<30(Xl,X2j = j ; i>i(Xl,X2j = 5^ i 02(Xl,X2j = 4 . (4) 

The Bispectrum B{k%, k 2 , kz) is defined in the Fourier domain as a three-point correlation function of the Fourier coefficients: 

<5(z)(ki)5(z)(k 2 )<5(.z)(k3)}c = (27r) 3 fc(k 1 + k 2 + k 3 )B s (fti, fe, fe, z). (5) 
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The skewness parameters are calculated by integrating over the bispectrum using the appropriate weights 
1 1 f°° kfdki f 1 



1 1 f°° k dki f 

S (k 2 ,z) = -r— -^r^T d f iB s (k 1 ,k 2 ,\k 1 +k 2 \,z)W{k 1 R)W(\k 1 + k 2 \R), (6) 
°o 4yr Jo 4yr J-i 

o -i poo 7,2 ji, pi 

Si(k 2 ,z) = - g 2 - . / / dM|ki+k 2 | 2 B i (fci,fc 2 ,|ki+k 2 | ! 2)W(fc lJ R)Vy(jk 1 + k 2 |7i), (7) 



9 1 rkfdkx f 1 , |2; 



ga(fel,g) = 4a 4 47r 2 _/ / d M (k 1 -k 2 )|k 1 +k 2 | / B 5 (fc 1 ,fc 2 ,|k 1 +k 2 |, Z ) 1 iy(fcli?)Ty(|k 1 +k 2 |i?), (8) 

/"OO 

S'i = / k 2 dkS\k)W(kR)- i =£ {0, 1, 2}. (9) 



with W^fci?) an arbitrary window function of smoothing radius _R. Although the window may be generic the two most commonly adopted forms 
are the top-hat and Gaussian window functions. We adopt the following notation: |ki + k 2 | = (k'f + k 2 + 2fcifc 2 /i) 1 / 2 and the angular terms are 
given by fi — (ki • k 2 )/(fcifc 2 ). In our convention the one-point skewness parameters can be recovered by integrating over the second momenta 
k 2 . For example: 

S (i) (z) = / k 2 2 dk 2 S (i) (k,z). (10) 



f°° 2 

(z) = / k 2 dk 2 %)(fc 
Jo 



In the approach presented above we are able to define a power spectrum associated to the MFs by considering the leading order corrections, 
O(ao), to the Gaussian MFs generalised to our skewness parameters: 



v { £ G) (v, k, z) oc {is (fc, z)H m+2 {u) + ^Si(k, z)H m {v) + m(m g 1] H m „ 2 (»)S 2 (k, z) J . (11) 

Essentially we have associated to each Minkowski Functional a power spectrum defined in terms of the three skew-spectra. The advantage 
to this approach is that we can study the contributions to the MFs as a function of the Fourier mode k. This is useful as a generic model for non- 
Gaussianity carries momenta dependence giving us greater distinguishing power than the conventional approach which compresses information 
from all modes into a single statistic. This compression represents a loss of distinguishing power. The skew-spectra presented above are relatively 
independent of our choice of background cosmology but will be dependent on the model for non-Gaussianity due to the aforementioned momenta 
dependence. 

Next we show how to generalize the skewness parameters to corresponding skew-spectra. 



4 NON-GAUSSIANITY: PRIMORDIAL AND GRAVITY INDUCED 

In the weak l y non -linear regime (8 < 1) the description of gravitational clustering can be described by Eulerian Perturbation Theory (see e.g. 
iBernardeaul 120021) ) . As the density contrast at a given scale becomes highly non-linear (5 > 1) the perturbative treatment breaks down and 
we observe a n increa s e in the growth of clustering. Perturbative studies of gravitational clustering have previously garnered a lot of attention. 
Starting with lPeeblesI Jl98Ct) there have been a significant number of attempts to reproduce the observed clusteri ng of a self-gravitating fluid in 
a cosmological setting, most of which adopt a brute force approach using N-body simulations (Bernardeau 2002). In the perturbative approach, 
solutions are generated by performing a series expansion with higher order corrections being introduced to the Fourier expansion of the linear 
density contrast under the assumption that the density contrast is less than unity for the series to be convergent: 

5(k) = «5 (1) (k) + «5 (2) (k) + «5 (3) (k) + ... ; 5 (2, (fc) = y^y^fo(ki+k a -k)F 2 (fci,fa)* ll) (ki)«W(k3). (12) 

The linearized solution for the density field is 5^ (k) with the higher order terms describing corrections to the linear term. Using a fluid approach, 
known to be valid at large scales, before shell crossing, one can write the second order corrections to the linearized density field by introducing 
a coupling kernel, F 2 (ki, k 2 ). Newtonian gravity coupled with the Euler and Continuity equations can be used to solve a system of non-linear 
coupled integral-differential equations, in order to generate the kernels F 2 (k\, k 2 ) and Fg(ki, k 2 , k^), by solving perturbatively order by order. 
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The expression for the matter bispectrum can be written in terms of an effective fitting formula that allows us to interpolate between the quasilinear 
and highly nonlinear regimes: 

B s (ki , k 2j k 3 ) = 2F 2 (ki , k 2 )P 5 (k x )P« (k 2 ) + cyc.perm. ; (13) 

2 



P 2 (ki, k 2 ) = ^-a(n e , ki)a(n e , fe 2 ) + ( k * , j*" 2 + k * , ^ 2 | 6(n e , fei)fe(n e , ft 2 ) + ^ ( -7— r-^ ) c ( n ei fci)c(n e , fc 2 ). (14) 
7 \ 2kf J 7 \ fci/c 2 / 

The coefficients a(n e , k), 6(n e , ft) and c(n e , k) are defined as follows: 

l + a 8 - - 2 ( 2 ),/( g /4)^ + 3.5 l + 0.4(n e +3)^+ 3 . f / 45_ 

a(ne ' fc) " 1 + ( 9 /4)«e+3.S > l + > C(« e ,fc)- 1 + (2?) „ e+3 . 6 |1+ ^ 15 + (ne+3)4 ^ 

(15) 

Here n e is the effective spectral slope associated with the linear power spectra n e = d In P a /d In ft, q is the ratio of a given length scale to the 
non-linear length scale q = k/k„t, where k 3 /2n 2 D 2 (z)Ps(k n l) = 1 and Q 3 (n e ) = (4 - 2" e )/(l + 2™ e ). Similarly cr 8 («) = P>(z)cr 8 . At 
scales where 5 < 1, and the relevant length scales are well within the quasilinear regime, then a — b — c = 1 and we recover the tree-level 
perturbative results. In the regime where 5 » 1, and the length scales under consideration are well within the nonlinear scales, we recover 
a = a^ ' 2 (z) \/0.7<33(n e ) with b — c = 0. In this limit the bispectrum becomes independent of configuration and we recover a hierarchical form 
for the bispectrum. The possibility, however, of weak violations of the hierarchical ansatz in the highly nonlinear regime is still not clear and can 
only be determined by higher resolution N-body simulations when they become available. Similar fitting functions for a dark energy dominated 
Universe calib rated against si mulations are also available and, at least in the quasilinear regime, most of the differences arise due to the linear 
growth factor dMa et alj|l999h . 

The analytical modeling of th e matter bispectrum pr esented here is equivalent to the the so called halo model predictions (see, for example, 
Mun shi. Smidt and Cooravl 1 2010h or lMunshi et al.l d20 1 lh for discussions of the perturbative treatment in the context of MFs). 



4.1 Primordial Non-Gaussianity: Bispectrum 

Much of the interest in primordial non-Gaussianity has focu sed on a phenomenological 'local /nl ' parametrization i n terms of the pertur bative 
non-linear coupling in the primordial curvature perturbation JSalopek and Bon dl 1 99d : iGangui 1994l: IVerde et al.ll2000l:lKomatsu et alJioOlh : 

= + /wi(*i(aO - <*£(a0» + 3NL$i(x) + h NL - 3<*i(x)» + . . . , (16) 

where <&l(x) denotes the linear Gaussian term of the Bardeen curvature and the amplitudes of the various non-Gaussian contributions are pa- 
rameterised by the set of variables {/nl, <?nl, ^nl, . . . }. In this parameterisation, the leading order non-Gaussian contributions are described by 
the bispectrum or, in configuration space, the three-point correlator. A large number of studies involvi ng primordial non-Gaussianity are based 
around the bispectr um as it will contain complete information (in a sta t istical sense) regarding /nl dBabichll2005l). The bispectrum has been 



around the bispectrum as it will contain complete information (in a statistical sense) regarding Jnl (riabich 2UU3). Ihe bispectrum has been 
extensively studied dKomatsu. Spergel & Wandelll2005l ; lc"rerninellil2003l : ICrermnelli et al j2 006: Med eiros & Co ntaldi 20061: ICabella et al. 1120061 : 



extensively studied ( Komatsu, Spergel & Wandelt 2UID; (Jreminelli 2 UUJI : HJrermnelli et al.l 2UUb; Med eiros & Co ntaldi 2UUa: Isabella et al. 2U0b; 
ISmith. Senatore & Zaldarriag3 l2009), with most of these measurements providing convolved estimates of the bispectrum. In Fourier space, the 
local-type bispectrum corresponding to[l6]has the following form: 

Bi oc (fci,fe 2 ,fei 2 ) = 2/^ L c [P*(fei)P*(fe) + cyc.] . (17) 

Standard inflationary models predict a primordial power spectrum that obeys a power law of the form: P*(fc) oc k n ~ 4 , where n is the spectral 
index. W e can explicitly relate the linear dens ity contrast to Bardeen's curvature perturbations on the matter dominated era through Poisson's 
equation (Hikage, Ko matsu & Ma stubara 20061): 

k 2 j> fc T(fc) = 47vG NPm (z) 7^7^ = ^n M H^S(k,z)(l + z). (18) 

In the above, T(k) is the tr ansfer function describing the ev olution of the density contrast. We adopt the approximation in lBardeen et al.l d 19861) 
(see Appendix B). See also lEisenstein and Hul d 19981 1 19991) for further discussions on transfer functions for large scale structure. At early times, 
where non-linear evolution may be neglected, the linear density contrast can be written in terms of an effective transfer kernel: 

50,, Z ) = ^=D( Z )M mk ; M( fc )4S4g. (19) 



(1+*) v-y-vv— . »vv 3n M H 2 ' 



© 0000 RAS, MNRAS 000, 000-000 



Non-Gaussianity in Large Scale Structure and Minkowski Functionals 7 



Here, D(z) is the linear growth factor normalised such that D(z) — > 1/(1 + z). This allows us to express the linear power spectrum in terms of 
the transfer kernel M(k) and the primordial power spectrum P<f,(k): 



Ps{k,z) 



M 2 (k) 



P*(k). 



(20) 



{i + zy 

We can relate the bispectrum for the primordial density perturbations to the bispectrum for the primordial potential perturbations by: 

Sr m (fci,fe 2 ,fei2,2) = D 3 (z)M(k 1 )M(k 2 )M(k 12 )Bl" m (k 1 ,k 2 ,k 12 ,z). (21) 

Therefore, in the linear regime (i.e. valid at large length scales), the primordial bispectrum for the local model B VTlm will evolve according to the 
following expression: 



B a ioc (k 1 k 1 ,k 3 ;z) = 2/^D 3 («) 



M(k 3 ) 



-P s (k 1 )P 5 (k 2 ) + eye. 



M(k{)M{k 2) 

The primordial potential bispectrum for the equilateral type can be expressed as JCremineW et al]2006h : 

B e q m = 6/ e r [_ ( p 4>(fcl)P$(fc2) + cyc .) _ 2(P 4> (k 1 )P $ (k 2 )P, ) (k 3 )) 2/3 + (P $ (k 1 )P|(k 2 )P|(k 3 ) + eye. 



1/3 



(22) 



(23) 



The primordial density bispectrum will therefore evolve according to: 



B, e9U! (k 1 k 1 ,k 3 ; 2 ) = 6/X^ 3 (« 



( M(k 3 ) 



VM(fci)M(fc 2 
-2 (M{ki)M(k 2 )M(k 3 )) 1/3 {P 5 (fc 1 )P i (fc 2 )P i (fc 3 )} 2/3 



)P s (k 1 )P s (k 2 ) +cyc.perm. 



M, 



1/3 



-{Psik^Psik^Psikaf}^ 3 + cyc.perm 



(24) 



Unlike the local model, the functional form for the equilateral model does not hold any relationship to fundamental physics but should be interpreted 
as something akin to a fitting function that describes a number of models for which the exact analytical expressions are more complicated. 
The folded configuration is maximised when two of the sides of the triangle are approximately equal such that k 2 « fc 3 w fcl/2. The folded 
bispectrum is approximated by the following functional form: 



B 



fold r> £ fold 

= 6/n L 



(P*(fci)P (fc 2 ) + cyc.) + 3(P< i ,(ki)P<i,(k 2 )P< 1> (k3)) 2/3 - (P$(ki)P|(k 2 )P|(k 3 ) + cyc.) 



1/3 



Following the outlined procedure we can re- write this in terms of the primordial density perturbations: 



B{ old (k 1 ,k 2 , k s ,z) = 



a ffold 
<> NL 

D(z) 



M(k 3 



-P s (k 1 ,z)P s (k 2 ,z) + cyc. 



M(ki)M(k 2 ) 

+ 3(M(k 1 )M{k 2 )M(k 3 ))- 1/3 {P s (k 1 ,z)P s (k 2 , z)P s (k 3 , z)} 2/3 
M^\k, 



M 1 /3(fe 2 )M(fc 3 



■{P s (k 1 ,z)P s {k 2 ,z) 2 P s (k 3l zf} 1/3 + cyc. 



(25) 



(26) 



The evolution of the primordial bispectrum is different compared to the one generated by gravitational evolution. The primordial bispectrum 
demonstrate momenta dependence and as a result the shape of the bispectrum configuration differs between the various primordial models as well 
as that for gravitational instability. At large angular scales, which will be probed by future weak lensing surveys, the gravitational instability may 
not have erased the memory of primordial non-G aussianity, which can provide complementary information to results obtained from CMB surveys 
such as Planck dThe Planck Collaborationll2006l) . 
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5 THREE-DIMENSIONAL DENSITY FIELDS 
5.1 Gravitationally Induced Bispectrum 

The first scenario that we will investigate is a three-dimensional smoothed field with galaxy density contrast S g . The matter bispectrum is given by 
Eq.lll3t which, when substituted into the skewness parameters Eq.® in the limit ijCl, gives: 



So{k2, z) 



kfdki 



dnW(kiR)W(k 3 R) 



'15 




r ki 


. k 2 \ 6 2 


_ 7 




\k 2 





P S {ki)Ps{k 2 ) 



Ol (2^) J (27T2) 

We now assume a Gaussian window function W{x) = exp(— x 2 /2) allowing us to make the following simplification (Matsubara 20021) : 

W (v + + ihhn) 



W{h)W{h) 



exp(—hl 2 /i); I = kR, 



(27) 



(28) 



with R the smoothing scale and k the Fourier mode. 

The angular integration over jj, can now be analytically performed by relating the terms in the above skewness parameter to terms corre- 
sponding to the m th Legendre p olynomial. This allows us to use the following relation between integrals over the Legendre polynomials and the 
modified Bessel functions J M (z) dLokas et all 1 9951 ; iMa tsubara 20021): 



m 2n 

-l)"y —Im + l/2(z). 



The relevant Legendre polynomials are 

P (z) = l; Pi(z) = z 
The resulting skewness parameters are: 

So{k 2 , z) 



P 2 {z) = \(Zz 2 -1); ft(z) = i(5* 3 -3*); 



Pi(z) = ^(35/ - 30z 2 + 3). 
8 



lldh 



Si(k 2y z) = 



S 2 (k 2 ,z) = 



sR^yf ww w2{h)wih)P6 (I) p * (I) ^ 



41 (h, h\ T nn ( 99 3 



a 4 (27r2)7 (27TW 



;2 ; 2' 

_L _|_ JL 
;2 + ;2 
'2 n . 

Zi 



W(Ii)W(fa)J% P<5 -5 V27r(ZiZ 



/3/a(iifa) + Ya (t 2 + |) - ^/r/ 2 (Ii/ 2 : 



\3/2 



^/ 1/2 (Z 1 Z 2 ) --^- + -J '3/2(^2) - -/ 5/2 (Z 1 Z 2 ) + ^ [- 2 + Ti ) Irrfhh) ~ ^(hh) 
where we have used the relationship Z = kR as introduced in Eq.l|28t. 

Using the notation for the skewness parameters proposed by Matsubara we can simplify the above equations with following variable: 



(29) 



(30) 



(31) 



(32) 



(33) 



S%f(fo, R) — V 



(7q 2tt 2 \aiR 



<T() 



a+0-2 



j^gimi 



(^-1^/2^01-3/2^-3/2 j 



i/2(Zi h 



(34) 
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3D Gravity Induced Non-Gaussianity in the Density Field z=0,1 




Figure 1. The three gravity induced skew-spectra So{k, z) (left panel), S\{k, z) (middle panel) and 52 (fc, z) (right panel) are depicted 3D cosmological density 
field. The spectra have been calculated for three different Gaussian smoothing scales R = 10, 20 h — *Mpc for z = 0. The skew-spectra for the redshifts z = and 
z = 1 are identical in the quasilinear regime due to the specific normalisation adopted. A ACDM background cosmology was assumed. The results are plotted for 
redshift z = 0. The skew-spectra are normalized in such a way that they are virtually independent of redshift. The expressions for So (k, z), S\ (k, z) and S2 (k, z) 
are given in Eq.gB, Eq.(32j and Eq.{33>. 



In this notation the skew-spectra parameters become: 

50 = y^ 1 -3S 1 02) + ^ 2 11 , (35) 

51 = S s « 13) - ^ 04) - lh 22 + t ^ 13) " S sf ' (36) 

a _ 81 c 33 9 (15) 99 33 9 (15) 72 33 

52 ~ 35 So " 5 Sl ~ 49* + 5 Ss " 245 S4 • (37) 

where we have adopted the following convention: S™^' = \{S?n + S^"). The results shown here demonstrate the analytic dependence of 
the gravitationall y induced b i spectr a on a weighted integral over the modes in terms of Bessel functions. These equations are analogous to 
those presented in iMatsubaral J2002l) but generalised to the skew-spectra formalism. The skew-spectra have been numerically calculated for three 
different smoothing scales i?[/i _1 Mpc] £ {10, 15, 20} and are shown in FigureQ] The non-Gaussianity generated from gravitational instability 
gives rise to a positively skewed per turbation unlike a generic inflationary model which could predict negatively skewed perturbations (e.g. 
Iffikage. Komatsu & Mastubara (2006)). It can also be seen that as the smoothing scale increases the amplitude increase but the scale at which we 
observe a cut-off incr eases (i.e. a cut-off occurs at smaller momenta). In the full skewness the gravitationally induced contributions are independent 
of smoothing scale l lHikage. Komatsu & Mastubarall 20061) . Additionally, due to the dependence of the gravitationally induced skew-spectra on the 
matter power spectrum and the variance of the field, we find that the generalised skew parameters are independent of redshift under the adopted 
normalisation. 



5.2 Primordial Bispectrum 

Using the above formalism we can calculate the skew-spectra for the primordial bispectrum in the local, equilateral and folded models. This is 
done by substituting the appropriate primordial bispectrum rEqs. J17l24l26H into the expressions for the skew-spectra rEqs. J6l7l8l >1 and evaluating. 
Due to the dependence of the primordial bispectra on the angular terms and the transfer function there does not appear to be a trivial analytical 
solution to these models and it is more elucidating to proceed by numerically integrating the resulting expressions allowing us to obtain the shape 
and amplitude of the spectra as a function of momenta. The resulting skew-spectra are shown in Figures ( 121314b . Unlike gravitational instability 
the primordial models have redshift dependence and exhibit similar behaviour to the gravitationall y induced models for varying smooth ing scales. 
In the full skewness the primordial models grow in amplitude as we increase the smoothing scale (Hikage, Komatsu & Mastubara 2006). 
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Local Primordial Non-Gaussianity in the Density Field z=0,1 
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Figure 2. The skew-spectra for primordial non-Gaussianity in a 3D cosmological density field is displayed for local type primordial non-Gaussianity. The spectra 
have been calculated for two different smoothing scales R = 10, 20h~ 1 Mpc with the left panel showing So(k, z), the central panel Si(k, z) and the right panel 
52 (fc, z). For each smoothing scale the upper curves correspond to higher redshift i.e. z = 1 and the lower curves correspond to z = 0. The local type primordial 
non-Gaussianity is considered which is defined in Eq. U7t . We have taken a ACDM cosmology with /ml = 20. 



Folded Primordial Non-Gaussianity in the Density Field z=0,1 
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Figure 3. Same as previous figure but for primordial non-Gaussianity of folded type Eq. d26t . 



6 THREE-DIMENSIONAL VELOCITY FIELD 

The linear density perturbations 8 g was introduced in Eq.dl7t as being related to Bardeen's curvature potential in the matter dominated era $. In the 
Newtonian limit the Universe can be treated as a pressure-less fluid and the continuity equation allows us to relate the linear density perturbations 
to an associated peculiar velocity v(k, z): 

4(k)+V-v(k) = 0. (38) 

In the basic inflationary paradigm and observationally the linear primordial density perturbations are highly Gaussian. Throug h the continuity 
equation we expect that deviations from Gaussianity will impact the observed peculiar velocity field. Following Matsubara (2002) we introduce a 
dimensionless scalar field normalised by the Hubble variable H: 

6(x) = H _1 V-v(x). (39) 
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Equilateral Primordial Non-Gaussianity in the Density Field z=0,1 
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Figure 4. Same as previous figure but for primordial non-Gaussianity of equilateral type Eo. i24t 

Gravity Induced Non-Gaussianity in the Velocity Divergence z= 0,1 
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Figure 5. The skew-spectra for gravitational instability in a 3D cosmological density field 0. The spectra have been calculated for three different smoothing scales 
R = 10, 20h Mpc with the left panel showing So(k, z), the central panel Si(k, z) and the right panel S2(k, z). The results for two different redshift slices are 
shown. For each smoothing scale the upper curves correspond to higher redshift i.e. z s = 1 and the lower curves correspond to z s = 0. The skew spectra are defined 
in Eq. )43L Eq. )44t and Eq.j45j respectively. The velocity divergence © is defined in Eq. )39t . The biaspectrum for is given in Eg. 1411 . The redshift dependence 
for velocity divergence bispectrum depends of gg defined in Eq. )42t . 



Due to the expansion of the Universe, the velocity field is described purely by the diverge (i.e. rotatio n -free) term with the rotational components 
being described by decaying solutions in perturbation theory (see, for example, Matsubara (2002); Bernardeau (2002)). Using second order 
perturbation theory the power spectrum and bispectrum can be analytically calculated dBernardeual dl994l) : 

P e (k)=g 2 g P s (k) + 0(at) 
"6 



Matsubara (2002)): 



Be(ki,k 2 ,k 3 ) 



-9 3 e 



, fci ,k 2 \ 8 2 

t 1 n^ + fct + r 



Ps(k 1 )P s (k 2 ) + cyc.perm. + 0(ctq). 



The factor go is the logarithmic derivative of the growth factor D(z) with respect to the scale factor a: 



dlnP 
d In a 



4/7 



A / ST2 
70 1 + -2 



M 



(40) 
(41) 

(42) 
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Figure 6. Same as previous figure but for non-Gaussianity induced due to primordial non-Gaussianity of local type with /jvl = 20. 




Figure 7. Same as previous figure but for non-Gaussianity induced due to primordial non-Gaussianity of equilateral type with /jvx = 20. 
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Folded Non-Gaussianity in the Velocity Divergence z=0,1 
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Figure 8. Same as previous figure but for non-Gaussianity induced due to primordial non-Gaussianity of folded type with /jvx = 20. 



6.1 Gravitationally Induced Bispectrum 

Using the modified power spectrum and bispectrum Eq. d41b we can calculate the skewness parameters for the peculiar velocity field: 
-1 1 f lldh 



So(h, z) 



Si(h,z) = 



S2H2, z) 



ge^o 



2n 2 R 3 



W\h)W{h)Ps 



geo-go-f An- 

33 fh 
28 U 2 



Ps 



X 



+ 



+ 



93 
35 



h/2{hh) + 



15 
14 



1 



ggaf 2tt 2 



I ) V2n(hl 2 ) 3/2 



12 



+ -)h/2(hh)-—l7/2(hh) 



(43) 



(44) 



WW2) 



9 



|v 2 (W 2 ) + A(| + | )/72!/|/j) 



144 
245 



-^9/2(^2) 



Using the parameterisation presented for the 3-dimensional cosmological density fields we can simplify the above expressions: 



(45) 



S (k 2 ,R) 



13 C n 3 (02) 



8 r,ll 
7 S2 



q (U T)\ - 1 ( 33 c.(!3) , 3 c,(04) 93 22 15 „(13) 12 22 

Sl (k 2 ,R)--(~-S +-S l --S, --S 2 +-5 3 



S2(fe)J? ) = i(_| S 33 + £ 5 r + |s| 3 . 

ge V 35 5 49 5 



9^(15) 



144 „ 33 
+ 245^ 



(46) 
(47) 
(48) 



These are, again, analogous to the appropriate expressions in iMatsubaral J2002h . The skew-spectra for the three different smoothing scales are 
shown in figure [5] 
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Gravity Induced Non-Gaussianity in Projected Density Field 




l/R (h Mpc ') 1/R (h Mpc ') l/K (h Mpc ! ) 



Figure 9. The skew-spectra induced by gravitational instability for projected (2D) cosmological density field. The spectra have been calculated for two different 
smoothing angular scales R = 10, 20/iMpc -1 with the left panel showing Co(l), the central panel Ci(/) and the right panel C2(l)- The solid line corresponds to 
R = 10/iMpc -1 and the dashed line corresponds to R = 20/iMpc -1 . The projected skew-spectra are defined in Eq. j56k Eq.^57) and Eq.^58). 



6.2 Primordial Bispectrum 

The 3D cosmological velocity field can be used to probe primordial contributions. We calculate the skew-spectra for the local model (Figure [6), 
the equilateral model (Figure[7} and the folded model (Figure[8j. All models have been calculated for redshifts z = 0, 1 and for smoothing scales 
R = 10, 20/i _1 Mpc. 



7 TWO-DIMENSIONAL PROJECTED DENSITY FIELDS 



Often cosmological data sets can be cast in the form of a projection over the sky. This requires us to introduc e a correspondenc e between three- 
dimensional density fields and the observed two-dimensional projection. The formalism adopted here follows Matsubara ( 2002). When working 
with a two-dimensional projected field one of the key observables are angular distances for objects placed at a given comoving distance r, Under 
the assumption of statistical homogeneity and isotropy the Universe can be approximated by the FRW cosmological models. The line element for 
an FRW cosmology can be expressed as: 



-cdt 2 + a(t) [dr 2 + d A (r) 2 (sin 2 ddd 2 + dcj> 2 )] 



(49) 



where we d A (r) is the angular diameter distance to radial comoving distance r. The scale factor is denoted by a(t). We have d A {r) 



sinh(K 



, r, sin(K 



for a Universe of negative, zero and positive curvature respectively. The curvature is given by K = (Qm — l)Hg. 



The three dimensional density field is then projected onto a two-sphere through the use of a selection function, n(r), and the angular diame- 
ter distance c?a(f) to relate the projected density contrast field to the underlying three-dimensional density field S g . Under the flat sky 
approximation this reduces to: 



dr d A (r) 5 g (r, r)n(r) 



(50) 



We can use a similar approach to move beyond one point statistics and consider the projected polyspectra. The two lowest order spectra of interest 
to us are the projected power spectrum P2D and the projected bispectrum Bid \ 



Pv{£) = J dr d 2 A (r) n 2 (r)P s 



f: 



d A (r)' 



By(l 1 ,£ 2 ,£ 3 )= / drd 2 A (r)n 3 (r)B 5 



£1 



d A (r) ' d A (r) ' d A (r) ' 



(51) 
(52) 
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Local Non-Gaussianity in Projected Density Field z=0,1 
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Figure 10. Same as previous figure but for local type primordial non-Gaussianity for /jvl = 200. 



Equilateral Non-Gaussianity in the Projected Density Field z=0,1 
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Figure 11. Same as previous figure but for equilateral type primordial non-Gaussianity. 



The three-dimensional power spectra Pg(k;r) has already been specified and the bispectrum Bg (ki ; r) will be evaluated for the appropriate model. 
In an analogous manner to the three-dimensional density fields we can construct variance and skewness parameters to allow us to investigate the 
contributions from various sources of non-Gaussianity to the projected bispectrum. The main difference in this approach is t hat we will be required 
to introduce the necessary formalism to cope with projection effects. The necessary machinery is discussed in iMatsubard h00$) though we have 
introduced the appropriate modifications to generalise the results to the proposed skew-spectra approach. The variance parameters of a projected 
field smoothed by an arbitrary window function, W(£8b), at a given smoothing scale, 9b, is given by: 



<Jj (0b 



2tt' 



eip{t)W 2 {ie b ) = 1 ^ [ drd 2 A (r)n 2 (r)D 2 (r)^[d A (r)e b ] 

"b 



We have introduced a variance function Ej to relate the three-dimensional field to our projected field: 

^(R) = R 2]+2 f ^k 2 >P l r{k)W 2 {kR). 



2tt 



(53) 



(54) 
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Figure 12. Same as previous figure but for folded type primordial non-Gaussianity. 



The skew-spectra parameters can then be found by introducing the appropriate functions of the projected field and projected derivative fields 
Eq.{§: 

S ln%) = {j^j 3 j drd 2 A {r)n\r)D\r)^- 2j [d A {r)6 b ) E? J [d A {r)6 b ]C^ [d A (r)8 b ] , (55) 

where skew-spectra parameters are given by: 

C (l 2 ,R) = ^ /^i / %=B 5 {h,h,h, z)W(h8 b )W(\h + h\6 b ), (56) 

Ci(l2,R) = =^ f 1 - 2 ^ 1 I %=\h+h\ 2 B s {h,h,h,z)W{he b )W{\\i + h\e b ), (57) 

^o^i J J 2n\/l — p 2 

C 2 (l2,R) = ^ f 1 -^ [ P -(li ■ l 2 )|li + h\ 2 B s (h, h, l 3 , z)W(hO b )W(\h + h\6 b ). (58) 

^1 J Ik J - fl 2 

For the 2-dimensional spectra the angular integration is performed by using the integral representation of the modified Bessel functions 
Im(hh) along with the recursion relations I' m = (J m -i + J m +i)/2. The integral representation is given by: 



i r 1 ili2M = fd_y = i f 1 d^ e „ ili2M 

W-i yrv V d{hh)j W-i yr^T? 

Proceeding analogously to the 3-dimensional case the 2-dimensional skew-spectra reduce to a weighted integral over the the modes k = l/R: 
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W 2 (h)W(l 2 ) 



93 

-jHhh) 



6 



3 r + r )h{hh) + -Hhh) 



7 



Ci(Z 2 ,fl) = 



2tt 



R 



p [-s) p [-s) w(h)W(l 



R 



§(i? + g)ib(W a )-| 



ll I 3 



27, 



C 2 (h,R) 



1 



hdl 



2tt 



fa + g + yiifa J Ji(JiIa) + 2g(i? + il)/ a (hi a ) 

2 



14 



hhhihh) 



W\h)W{h) 



^-llllh{hh) ~ ~{l\h + llli)Ii(hh) - Y l i%Hhh) - + llh)h(hh) + lillh{hh) 



7 x " v ' 2 
We adopt a parameterisation in terms of the variable C^: 

a + P-2 



C%f(h, R) — ^4 
^0 

The skew-spectra can then be written as: 



hdl 



2tt 



(60) 



(61) 



(62) 



(63) 



C (h,R)= 2 4c?-6Cr + 6 -C 1 2 \ 



7 

C 1 (l 2 ,R)= 9 ^ct ) 



7 

(40) 8 1 22 . 33 ((31) 3 22 

3Cl ~ 14 Gl + 14 62 ~ 14 63 ' 



C 2 (Z2, fl) = y C7 33 - 3C{ 42) - y C 2 33 - 3C? + Cf. 



(64) 
(65) 
(66) 



The resulting skew-spectra for gravitational instability are shown in Figure l[9} for smoothing scales R = 10, 20 frMpc -1 . 

The 2D projected surveys can be used with the primordial bispectra by simply replacing the bispectrum kernel in Eqs. d56l57l58t with the 
corresponding primordial model. The results are shown in Figures d 101 1 1112V 



8 REAL WORLD EFFECTS 

In this section we briefly outline some key observational and systematic effects that lead to corrections and modifications to measured quantities 
(e.g. power spectrum). This paper does not attempt to deal with real world effects but a more general treatment incorporating real world effects 
will be presented elsewhere. 



8.1 Redshift Space Distortions 

Real world observations and surveys consist of constructing a map detailing the frequency and intensity of radiation across the sky or a patch of the 
sky. Cosmological distances are constructed under the assumption that the frequency of a particular emitter (e.g. CO emission) is known allowing 
us to estimate the redshift to the source. Under the assumption that we inhabit a spatially homogeneous, isotropic, irrotational and expanding 
spacetime then the redshift may be converted to a coordinate distance in a fairly straightforward manner. Redshift distance estimates consist of a 
number of contributions of which the two most physically relevant, for cosmological observations, are the metric expansion of spacetime (as the 
dominant factor) and the Doppler shifts arising from local peculiar velocities of a source along the line of sight. 

The non-linear nature of the mapping from real space to redshift space means that a Gaussian random field will have non-Gaussian corrections 
when observed in redshift space. In addition to these non-Gaussian corrections there will be modifications to the observed power spectrum. Such 
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system atics will become increasingly importan t with the precision of the surveys mentioned in Section[T] There exists analytic results for the linear 
results jKaiseJ Jl987t);lHamiltonl ll992lll997l) : lMatsubard J200C|) 1 and recent studies have begun to investigate the full non-linear redshift-space 
power spectra (e.g. lShaw and Lew is (2008)). The linear result first provided in lKaiserl<ll987h introduces the following modification to the 2-point 
function: 

P g , t (k) = b\l + b^ffilfPsik), (67) 
where b is a linear bias factor, f the derivative of the linear growth factor and /j,k = n ■ k. 

8.2 Galaxy Bias 

The formation of galaxies and collapsed objects is highly dependent on the total local matter density. An important point to note is that large 
scale structure surveys only probe the luminous baryonic matter and we find that there are nonlinearities in the biased relationship between 
the observable baryonic matter distribution (i.e. galaxies) and the underlying total matter distribution (i.e. baryonic plus cold dark matter). This 
nonlinearit y can give rise to an effect ive non-Gaussianity that must be taken into account when performing statistical tests of large scale structure. 
Following iFrv and Gaztanag j fl993h we can assume that galaxy formation is a local deterministic process and perturbatively expand the galaxy 
density contrast S g in terms of the underlying matter density contrast 8 m : 

8 g {z) = b (z) + b 1 (z)S m (z) + ^S 2 m (z) + O(Sl), (68) 

where bo(z) is set by the constraint (S g (z)) — 0. The two remaining parameters b\(z) and 62(2) are the galaxy bias parameters. The 2-point 
and 3-point correlators of the galaxy distribution are therefore weighted functions of the 2-point and 3-point correlators of the matter distribution: 

P g {k,z) = bl{z)P m {k,z), (69) 
B g (ki,k 2 ,k3,z) = b\{z)B m {ki,k 2 ,k 3 ,z) + bl{z)b 2 {z) [P m (ki, z)P m {k 2 , z) + eye] . (70) 

In the above B m (ki, k2,k^, z) is the bispectrum of the matter bispectrum. For this report we will neglect the bispectral contribu t ions o f 
galaxy biasing and focus prima r ily on the pr i mordial non-Gaussianities, for a detailed discussion of biasi ng effects see lFrv and Gaztanagal £l993); 
IScoccimarol j2()00): McDonald 1 20oA l2008h : ISefusatti and Komatsul l2007h : |jeong and Komatsul [2009iM for a selection of relevant papers. 

8.3 Survey Effects 

In addition to redshift space distortions and biasing effects we also have to consider additional observational effects such as survey geometry, survey 
volume and corrections to the radial selection function arising from magnitude limited surveys. Other effects include the survey specifications (e.g. 
range of redshifts probed), instrumental noise and additional complex non-linearities that can manifest at high k altering the power spectrum. 
Examples of these complex non-linearities include both instrumental transfer function effects in Fourier space and survey specific sample biasing. 
A specific model for the instrument transfer function must be assumed to investigate effects at high k for a given survey. Such a modification is 
arbitrary and survey dependent. Skew spectra, being high order statistics, will be affected more by sample (cosmic) variance. The resulting bias 
and scatter of these estimators will have to be dealt with more precautions than their ordinary power spectrum counterparts. Such issues will be 
dealt with in a separate publication. We shall provide a quick dis cussion of magnitude limited surveys, radial selection function and volume of the 
survey following the discussion presented in (Hikage et al. 20q!). 

In an apparent-magnitude limited survey the number density of observed galaxies will decrease with redshift. This means that for the galaxies 
to be appropriately resolved, and hence a reliable density field reconstructed, we need an appropriate smoothing length R > R g such that R g 
is the mean separation of galaxies at our specified redshift. This means that we have a maxim um redshift z max = Zma, x (R g ) such that R g will 
be the mean separation of the magnitude limited galaxy sa mple at our maximu m redshift z max dVogelev et al 1 19941) . The number of independent 
resolution elements in each survey volume will be given by 1 Hik age et alj2Q02h A r R OS = V S urv(< Z max )/(2n) 3 ^ 2 R g . This paper does not attempt 
to take real-world systematics into account, a treatment of real-world systematics and error estimates will be presented elsewhere. 
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9 CONCLUSION 

New surveys proposed over the next decade or two offer the possibility of probing large scale structure in an unprecedented manner. It is hoped 
that larger, more comprehensive cosmological data sets along with a more sophisticated understanding of the underlying systematics will allow 
us to investigate the role of non-Gaussian initial conditions along with the physics underlying gravitational instability in an FLRW Universe. This 
makes the statistics of galaxy clustering an interesting playground for model testing and discrimination. 

In this paper we have adopted the conventional approach to describing gravitational clustering through a hierarchy of higher order correlation 
functions. There exist numerous statistical methods to characterise the higher order correlators ranging from polyspectra defined in Fourier space 
to the spatially defined Minkowski Functionals adopted in this paper. The MFs are a topological statistic that characterise the morphological 
properties of a random field. In particular it has been shown that the morphological properties of a Gaussian random field are analytically known 
with non-Gaussian corrections giving rise to changes in the morphological properties of a random field. The MFs are therefore sensitive to non- 
Gaussianity allowing us t o cons truct a statistical estimator based on a nalyti cally known results for a Gaussian random field. Using the formalism 
presented by Matsu baral J2002h and iHikage. Komatsu & Mastubaral d2006h it has been shown that the MFs, to leading order, depend on a set 
of skewness parameters So, Si, S2 that encapsulate the non-Gaussian corrections to the random field. In particular, for a weakly non-Gaussian 
random field, the skewness parameters, a set of one-point estimators, reduce to a weighted probe of the bispectrum. The currently adopted approach 
has been to collapse all the information into the one-point estimators but at a loss of the ability to effectively discriminate between different types 
of non-Gaussianity. The shape of the bispectrum has been strongly emphasized as an extremely important tool in discriminating between various 
models for an inflationary scenario as well as secondary non-Gaussianity such as gravitationally induced non-linearity. Motivated by this we have 
presented a generalisation of these one-point estimators to a set of skew-spectra Sf , Sf , Sf that avoid collapsing all the information into a single 
estimator. It has been argued that by investigating the shape of the skew-spectra and it's relation to the bispectrum Bi t i 2 i 3 we can discriminate 
more effectively between the various bispectral configurations and ideally separate out the primordial and secondary contributions. 

These generalised skew-spectra are two-point statistics spatially but the dependence on the bispectrum is related to the fact that they are third 
(leading terms) order statistics in terms of non-Gaussianity. A cumulant correlator of order p + q is constructed from the cro ss-correlation of the 
spatia l fields (or the derivative fields) (4 ,p (x)^I' 9 (x)). This cumulant correlator will probe the polyspectra of order p + q jMunshi & Heavens! 
l2009h . 

This paper has only considered the generalised skewness parameters up to O(ao) with the resulting spectra corresponding to a weighted 
probe of bispectrum. The formalism is developed in terms of a perturbative expansion about a Gaussian random field and we necessarily truncate 
the results to the desired order. We could consider highe r order corr e ctions with the resu lts of 0((7q) for which the resulting spectra would 
provide a weighted probe of the trispectrum (see e.g. lMatsubaral ( 2010); M unshi et al.l J201 lh ). The next-to-leading order results are appropriately 
named kurt-spectra due to their relation to the more conventionally used kurtosis parameters. Although we have not considered the next-to-leading 
order non-Gaussian corrections, the unprecedented level of detail in upcoming surveys makes investigations of next-to-leading order statistics a 
promising future direction. 

The skew-spectra have been constructed from various products of the spatial and derivative fields (e.g. V>t, V 2 \l/, • V\l/, etc). The 
skew-spectra can be related to the spatially defined cumulant correlators and the skewness parameters can be constructed from the corresponding 
skew-spectra. Given the momenta-dependence of the skew-spectra it was shown that they carry greater discriminating power against the various 
bispectra configurations which can be compared to observational or numerical results. Something that has not yet been considered is the role of 
noisy data, systematics or survey masks. A discussion of error on the estimators and their scatter will be presented elsewhere. The pre sence of such 
noise will lead to scatter in our estimators necessitating a greater understanding of the systematics in upcoming LSS surveys (see e.g. lllikage et al.l 

The systematic study of the skew-spectra presented in this paper focused on three different cosmological data sets: 3D density, 3D velocity 
and 2D projected. For each of these models we numerically investigated the shape of the bispectrum for three primordial models (local, equilateral 
and folded) along with the analytical gravitationally induced bispectrum. The primordial bispectra were shown to be dependent on the underlying 
matter power spectrum, transfer function with an overall amplitude set by the background cosmology and smoothing scale. The 2D skew-spectra 
were calculated for the innermost kernel in the integral with the integral over the background geometry neglected. The shape and amplitude of the 
underlying spectra is still of interest however. 

The approach to probing non-Gaussianity presented in this paper is just one of many areas of similar work. It is accepted that the CMB 
should provide one of the cleanest probes of the early Universe with the density perturbations being adequately modeled by linear perturbation 
theory. The non-linear nature of gravitational instability means that the primary source of non-Gaussianity in LSS or weak-lensing surveys is most 
likely non-primordial. Nevertheless, LSS surveys are currently thought to be able to place comparable constraints on primordial non-Gaussianity 
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to that of the CMB. Likewise, recent developments on the analytical modeling of weak lensing has lead to the hope that upcoming weak lensing 
surveys could provide an unbiased pr obe of the underlyin g matter distribution yielding a cleaner probe of the statistics of gravitational clustering 
than conventional LSS surveys (e.g. lMunshietalJ < l201ll) ). Another possibility that could offer the cleanest probe of primordial non-Gaussianity 
is that of a non-Gaussian stochastic background of gravitat ional waves. The pros pects for direct detection are not currently optimistic due to the 
weakness of the 3-point term in the graviton interactions dAdshead and Limll2009h though it is hoped that indirect detection through, for example, 
CMB polarisation could be a possibility. 

An assumption that was made throughout this paper was the use of a Gaussian window function. This is a popular, though arbitrary, choice 
and we could equally have p erformed the analy sis using alternative window functions. Another popular window function used in the literature is 
the top-hat window function dBernarde au 2002). The formalism presented allows the results to be generalised to arbitrary functions but this is not 
thought to change the overall conclusions of the paper. 

In all models the gravitational instability bispectrum was an order of magnitude or two greater than that of the primordial bispectra. At 
very small wavenumbers the local primordial bispectra was starting to become comparable to the gravitational instability skewness. The results 
quoted are for /nl ~ 20 which may be somewhat optimistic with the true amplitude for non-Gaussianity lying closer to unity and as a result this 
would further reduce the skew-spectra amplitudes by an order of magnitude. Interestingly it seems that the shapes of the primordial bispectra are 
distinguishable even if the gravitational instability non-Gaussianity is the dominant contribution for all surveys discussed. The results do not take 
into account systematics or noise so no constraints on the signal-to-noise ratio may be discussed but a more detailed account will be presented 
elsewhere. The shape of the gravitational instability skew-spectra could be useful in quantifying predictions from perturbative treatments of 
gravitational instability, N-body simulations or semi-analytic approximations and it is worthwhile considering the prospects for disentangling any 
primordial contributions from gravitational instability. 

We have studied the skew-spectra for the 3D galaxy distribution as well as 2D projected surveys. For the 3D surveys we ignored the ef- 
fects of redshift space distortion. Redshift space effects are important as these correspond more directly to the observables in conventional large 
scale structure surveys. Examples of the role of red shift space distortion can be found in jKaiseJI 19871 : lHamiltorj[l997l : | Shaw and Lewisll2008l ; 
Oku mara and Jin j|201 it [Sato and Matsubara 201 2). A more general treatment taking into account redshift space effects will be presented else- 
where. 
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APPENDIX A: SPECIAL FUNCTIONS AND THEIR PROPERTIES 

A number of calculations in both this paper and key references make use of a number of properties of the special functions (e.g. the Legendre 
polynomials). In this section we provide some key results used in the derivations. 

In order to perform the angular integration in the skewness parameters it is convenient to use the relationship between the Legendre poly- 
nomials and the modified Bessel fun ctions. The origin of this relationship can be seen by performing an expansion of the angularly dependent 
Gaussian window function W{kzR) dLokas et al.ll 19951) : 



g-li-l2 _ g-Vap 



Y, M)" 1 (2m + l) A /^/ m+1/2 (W 2 )P m (M)- (AD 



The Legendre polynomials P m (/i) inherit the angular dependence and upon integration we can use the orthogonality and completeness of these 
functions to simplify our results considerably: 

1 ^P m (^)P n ( M ) = * &m 2,y (A2) 

.1 (2m + 1) 
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This relationship motivates us to associate the angular term s dependent on fi to the appropriate Legendre polynomial. After a small amount of 
algebra we recover the formula provided in iMatsubaral f2002h relating the Legendre polynomials and modified Bessel functions: 

/ dnP m (n)e- hl ^ = (-1)» A /Ill m+1/2 (W 2 ). (A3) 
J-i V nfo 

A useful set of relationships are those for the 2-dimens ional case in whic h the recursion relationships allow us to perform the angular integration. 
The following formulas are particularly useful (see e.g. lMats ubara (20Q2|)): 

1 I "'' -_e~ 1 ^ =I (hh), (A4) 





















/< 2 


d/J, 














/< 2 



/xe-' 1 '^ = -Ii(Ii i a ), (A5) 



- [ -J= — _ /i 2 e~' 1 ' 2A ' = i (Jo(«i/ 2 ) + / 2 (l!l 2 )) , (A6) 

71" ./ a /I — II 2 2 



' "" = M 4 e~ ili2A1 = i(3/oai/2)+4/ 2 (/ 1 Z 2 ) + / 4 ai/2)). (A8) 

These results were used in the text of the paper. 

APPENDIX B: TRANSFER FUNCTION 

The transfer function allows us to evolve an initial sub-Hubble mode to some final state at an arbitrary time. 

T ^ = ^T7TrT (B1) 

o(k,ai) D(a) 

The growth function defines the linear growth of long wavelength perturbations ( lBardeenetaT]|l986h . We require that the initial time is such that 
it is before any scale of interest has become sub-Hubble. By definition we have: 

$(fc,a) =T(k,a)<j)(k,a i ) (B2) 

This allows us to evolve the initially sub-Hubble modes from a given model for the early-Universe to make predictions for the power spectrum 
at an arbitrary time given that the approximation for the transfer function holds. The model we adopt for the transfer function is that given in 
teardeenetatlll986h for a cold dark matter model with adiabatic initial conditions: 

T(q) = lnl 2 + 3 y 4g [1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71#] ~ 1/4 (B3) 
where we have adopted the following variables: 

j-; T = a m , Q h 2 (B4) 



(rMpc- 1 ) 

This fitting function built on the work of ( Bard ee3l985t iBond and SzaTavll 19831: | Bond and Efstafhioull 19841 : Efstathiou and Bonj|l985t) . 
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